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We report measurements of oscillation parameters from i/ M and 77 M disappearance using beam 
and atmospheric data from MINOS. The data comprise exposures of 10.71 x 10 20 protons on tar- 
get (POT) in the ^-dominated beam, 3.36 x 10 20 POT in the ^-enhanced beam, and 37.88 
kton- years of atmospheric neutrinos. Assuming identical v and 77 oscillation parameters, we 
measure | Am 2 1 = (2.41^;^) x 10" 3 eV 2 and sin 2 (2fl) = 0.950lg;o36- Allowing independent v 
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and v oscillations, we measure antineutrino parameters of |Am | = (2.50_ ' 2 5 ) x Hn^eV and 
sin 2 (26) = 0.97io Qg, with minimal change to the neutrino parameters. 

PACS numbers: 14.60.Pq, 14.60.Lm, 29.27.-a 



Neutrino oscillation provides direct evidence that neu- 
trinos have non-zero mass and represents the only phe- 
nomenon observed to date with an origin beyond the 
Standard Model of particle interactions. With massive 
neutrinos, three flavor eigenstates mix with three mass 
cigenstatcs according to a unitary matrix that can be pa- 
rameterized by three angles and a CP-violating phase [lj . 
The resulting oscillation probability depends on the mix- 
ing angles and on the differences between the squared 
neutrino masses. The MINOS experiment performs pre- 
cision measurements of oscillations via disappearance. 
These oscillations are well described by an effective two- 
flavor model with flavor and mass eigenstates related by 
a single mixing angle 6. In this approximation, the 
survival probability is given by 

, , . 2 „„ . , / 1.267 Am 2 [eV 2 ] L[km 

P iy u v u ) = l-sin 2 20 sm 2 — ^ — l - ■ 

\ E[GcV] 

where L is the distance traveled by the neutrino and E is 
its energy. The mass splitting, Am 2 , is an admixture of 
the three- flavor parameters Am 2 ^ and Am| 2 0], and it 
differs from Ato 2 2 by less than 2%. The survival prob- 
ability has the same form, but the mixing parameters are 
denoted by Am 2 and sin 2 (26) . 

The MINOS measurements use neutrinos produced in 
the NuMI accelerator complex and by cosmic ray inter- 
actions in the atmosphere. The accelerator provides a 
source of neutrinos with a fixed baseline and an energy 
spectrum that peaks at L/E ~ 250km/GeV, close to the 
region where the survival probability reaches its first 
minimum. Atmospheric neutrinos are produced with a 
broad range of£7(~ 0.5-10 4 GcV) and L(~ 10-10 4 km), 
enabling the study of oscillations across a wide region in 
L/E and covering many oscillation cycles. The precision 
of the oscillation measurement is enhanced by combining 
these two complementary samples. 

This Letter presents the first ever joint analysis of at- 
mospheric and accelerator neutrinos in the same experi- 
ment. The new results come from the full MINOS data 
set, collected over a period of nine years. The combi- 
nation of data sets, together with increased exposures, 
produces a significant improvement in the sensitivity to 
oscillations over previous MINOS analyses Fur- 
thermore, MINOS has the unique ability to separate neu- 
trinos and antineutrinos on an event-by-event basis. Cou- 
pled with the world's only set of long-baseline accelerator 
antineutrino data, we present the most precise measure- 
ments to date of the larger mass splitting for both neu- 
trinos and antineutrinos. 

The NuMI beam Q is produced at the Fermi National 
Accelerator Laboratory (Fermilab) by 120 GeV protons 



striking a graphite target. The resulting charged pions 
and kaons are focused by two magnetic horns before de- 
caying in a 675 m long helium- filled volume 0] . The beam 
is directed through a hadron absorber and rock to stop all 
particles except neutrinos. The energy spectrum of the 
neutrino beam can be changed by varying the distance 
between the target and first horn. Most of the data used 
in this analysis were collected with a spectrum peaking 
at a neutrino energy of 3 GeV. By selectively focusing 
positive or negative pions and kaons, a ^-dominated or 
77^-enhanced beam is produced. 

The two MINOS detectors are steel-scintillator track- 
ing calorimeters with toroidal magnetic fields Each 
detector consists of steel plates with segmented plastic 
scintillator planes mounted on them. The planes are per- 
pendicular to the beam direction. The 0.98 kton Near 
Detector (ND), located 1.04km from the neutrino pro- 
duction target, measures the beam composition and en- 
ergy near the source. The 5.4 kton Far Detector (FD) 
measures the beam composition and energy spectrum 
735 km away from the target. Installed 705 m (2070 m 
water-equivalent) underground in the Soudan Under- 
ground Laboratory in Minnesota, the FD is also used 
to measure oscillations in atmospheric neutrinos and an- 
tineutrinos. A scintillator veto shield is installed above 
the FD to enhance the rejection of the cosmic-ray muon 
background. 

Muon neutrinos and antineutrinos arc identified 
through their charged current (CC) interactions 

The muon typically deposits energy in the detector in a 
clear track-like pattern. The hadronic recoil system, X', 
leaves a diffuse shower-like deposition pattern. The only 
notable background in the CC sample arises from a small 
number of neutral current (NC) interactions that gener- 
ate only hadronic activity but can display a track-like sig- 
nature. Muon neutrinos and antineutrinos are separated 
by the direction of curvature of the charged muon track 
in the magnetic field of the detectors. The muon momen- 
tum is determined from the range for muons that stop in 
the detector and from curvature for exiting muons. For 
beam neutrino interactions, a /c-Nearest-Neighbor classi- 
fication algorithm (fcNN) is used to estimate the hadronic 
energy from both the calorimetric energy deposited and 
the topology of the shower [lo| . For atmospheric neu- 
trino interactions, the calorimetric energy deposits in 
each scintillator strip are summed to provide an estimate 
of the true shower energy. The reconstructed neutrino 
energy is given by the sum of the muon and shower en- 
ergy measurements. 
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dominated beam and 3.36 x 10 20 POT in the F p -enhanced 
beam, corresponding to increases of 48% and 14%, re- 
spectively, over our previous analyses @, 0] . As in these 
previous analyses, the selection of and F M CC interac- 
tion candidates proceeds via the construction of a set of 
variables that characterize the event topology and energy 
deposition of muon tracks Again, these variables 



are combined into a single discriminating variable using 
a fcNN technique. 

From the ^-dominated beam, we use both neutri- 
nos and antineutrinos with interaction vertices contained 
within the detectors' fiducial volumes. When explicitly 
fitting antincutrino oscillation parameters, we apply an 
optimized antincutrino event selection to increase the pu- 
rity of this contained- vertex antineutrino sample Q. In 
this beam mode, we also select a sample of non-fiducial 
muons in the FD, comprising CC and CC inter- 
actions outside the fiducial volume and in the rock sur- 



rounding the detector [FJ, LL3[ . No muon charge-sign sep- 
aration is performed on this sample since many of these 
muons are confined to the edges of the detector where the 
magnetic field is very low and muon curvature is less well 
modeled. In the F^-enhanced beam, only the contained- 
vertex antincutrino sample is used, as the non-fiducial 
sample is dominated by high energy neutrinos. 

The predicted FD beam spectra are derived from the 
observed ND beam data using a beam transfer ma- 
trix [l4[. This extrapolation procedure mitigates many 
sources of systematic uncertainties and naturally ac- 
counts for any variations in the beam conditions such as 
target degradation or differences among the seven differ- 
ent production targets used throughout the experiment's 
lifetime. Since the ND is used to provide a baseline spec- 
trum, it is important to minimize any differences between 
the response in the two detectors. In particular, the re- 
gion around the ND magnetic coil is poorly modeled, so 
any beam-induced events with muon tracks entering this 
region are removed from the ND data set. 

We use the same atmospheric neutrino data set and 
event samples as our previous analysis [H , which is based 
on a FD exposure of 37.88 kton-ycars. The events are 
identified by the presence of either an interaction vertex 
within the fiducial volume of the detector or an upward- 
going or horizontal muon track. The selected events are 
separated into three samples: contained-vertex muons, 
non-fiducial muons, and contained-vertex showers. The 
two muon samples are produced by CC and CC 
interactions; the contained-vertex shower sample is com- 
posed mainly of v e CC, V e CC and NC interactions. The 
atmospheric neutrino samples must be selected from a 
background of cosmic-ray muons. For contained-vertex 
muons and showers, this background is reduced to 4% 
and 12%, respectively, by applying a series of contain- 
ment requirements and by checking for energy deposits 



Data Set 


No osc. 


With osc. 


Observed 


Vy, from beam 


3201 


2543 


2579 


V : M from Vf,, beam 


363 


324 


312 


Non-fiducial fj, from beam 


3197 


2862 


2911 


77 M from 77 M beam 


313 


227 


226 


Atm. contained-vertex i/ M + 77 M 


1100 


881 


905 


Atm. non-fiducial fi~ + 


570 


467 


466 


Atm. showers 


727 


724 


701 



TABLE I: Numbers of events selected in each sample. The 
oscillated event yields come from the best fit to all data, as- 
suming identical u and V oscillations. 



in the sections of veto shield above the fiducial event ver- 
tex within a ±50 ns window. In the non-fiducial sample, 
the background is almost entirely removed by using the 
2.5 ns timing resolution to accurately determine the in- 
coming muon track direction. TableUlists the numbers of 
observed events and the corresponding predictions, with 
and without oscillations, for each of the analyzed sam- 
ples. 

We simulate atmospheric neutrino events according to 
the Bartol flux calculations [l5| . The beam neutrino flux 
is simulated using the FLUGG [l6| package, which com- 
bines GEANT4 [nl geometry with FLUKA [HI hadron 
production. All beam neutrino interactions, and inter- 
actions of atmospheric neutrinos in the detectors, are 
simulated using NEUGEN3 [l8|. We simulate atmo- 
spheric neutrino interactions in the surrounding rock 
using NUANCE [l9[, which propagates the final-state 
muons through the rock and up to the edge of the de- 
tector. In both the beam and atmospheric simulations, 
the propagation of particles in the detector, and the de- 
tector response, are simulated with GCALOR [2(} and 
GEANT3 " 



211 ] . The simulation incorporates the back- 



ground arising from v ^ — > v T appearance. For the best 
fit oscillations, the predicted event yield from this chan- 
nel totals 18 events across the entire data set. 

The oscillation parameters are obtained from a maxi- 
mum likelihood fit to the data. The measured FD beam 
data are binned as a function of reconstructed neutrino 
energy To improve the sensitivity of the analysis, the 
contained-vertex events from the ^-dominated beam 
are divided into five sub-samples according to their es- 
timated energy resolution, which lies primarily in the 



range 5 — 30% [22j, l23j. The atmospheric samples are 
binned as a function of L/E. The contained-vertex at- 
mospheric and v ^ events are divided into four sub- 
samples according to the estimated L[km]/i5[GeV] res- 
olution, where 0i oglo (l/e) ranges from 0.05 to 1.2 (24|. 
The contained-vertex showers are grouped in a single bin 
because the majority are too short for an accurate mea- 
surement of L/E. These events are relatively insensitive 
to oscillations but provide a constraint on the overall nor- 
malization of the atmospheric flux. 
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FIG. 1: FD data samples compared to predictions with and without oscillations. The top row shows the energy spectra of the 
beam samples, while the bottom row shows the L/E distributions for the atmospheric event samples. 



The fit incorporates a set of nuisance parameters that 
accommodate the largest systematic uncertainties in the 
simulation of the beam 0, [H| and the atmospheric Q neu- 
trino data. For both data sets, the fit incorporates the 
systematic uncertainties in the overall normalizations of 
the event samples, the relative normalization of the NC 
background component, the muon momentum, and the 
shower energy. The latter two uncertainties are taken as 
correlated between the beam and atmospheric samples. 
An analysis performed with all uncertainties uncorrelated 
produces similar results. Additional systematic parame- 
ters are included in the fit to cover the uncertainties in 
the rate and spectral shape of atmospheric and 
events arising from uncertainties in the neutrino flux and 
cross-section simulations. 

When wc fit the full MINOS data sample to the two- 
flavor neutrino oscillation hypothesis, assuming that neu- 
trinos and antineutrinos have identical oscillation pa- 
rameters, we obtain |Am 2 | = (2.41+g;?^) x 10~ 3 eV 2 
and sin 2 (26>) = 0.950±o;o3g. Maximal mixing is disfa- 
vored at the 86% confidence level (C.L.); we measure 
sin 2 (26>) > 0.890 at 90% C.L. The observed beam and 
atmospheric event spectra in the FD are shown in Fig. [T] 
along with the predictions for the case of no oscillations 
and the best fit. The data are well described by the neu- 
trino oscillation model; the same analysis performed on 
simulated experiments returns a worse quality of fit for 



19.1% of those experiments. A number of cross checks 
were performed by fitting each of the data samples sep- 
arately. Those separate fits yielded consistent oscillation 
parameters, indicating that the data samples are consis- 
tent with each other and with the oscillation hypothesis. 
Allowed regions for the oscillation parameters, assum- 
ing identical neutrino and antincutrino oscillations, are 
shown in Fig. [2] 

The magnetized MINOS detectors enable separation 
of neutrino and antincutrino interactions for both 
beam and atmospheric events, allowing an independent 
measurement of the antincutrino oscillation parame- 
ters. We perform an additional fit in which wc allow 
neutrinos and antineutrinos to have different oscillation 
parameters, and find |Ato 2 | = (2.50±£fJ) x 10" 3 eV 2 
and sin 2 (26) = 0.97±g;g| (> 0.83 at 90% C.L.). 

The difference between the antineutrino and 
neutrino mass splittings is measured to be 
|Am 2 | - |Am 2 | = (0.12±g;||) x 10- 3 eV 2 . Correspond- 
ing measurements using the beam and atmospheric 
samples separately yield consistent results. The 
90% C.L. allowed region for the antineutrino oscillation 
parameters is shown in Fig. |3l illustrating good agree- 
ment between the measured neutrino and antincutrino 
oscillation parameters. 

In summary, wc have presented an analysis of the 
combined MINOS beam and atmospheric neutrino sam- 
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FIG. 2: The allowed regions of |Am 2 | and sin 2 (2#), assuming 
identical neutrino and antineutrino oscillations. The MINOS 
result is compared to results from Super-Kamiokande [25| and 
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FIG. 3: The 90% confidence level allowed region of |Am 2 | 
and sin 2 (28) from the fit assuming identical neutrino and an- 
tineutrino oscillations (MINOS 90%, also in Fig. [2| is com- 
pared to the allowed regions for and 77 M from the fit in 
which neutrinos and antineutrinos have different oscillation 
parameters. 



pies, representing the complete data set from the MI 
NOS experiment. Assuming that neutrinos and antineu 
trinos share identical oscillation parameters, we mea- 
sure sin 2 (26») = 0.950^;^ (> 0.890 at 90% C.L/ 
and |Am 2 | = (2.41±°;^) x 10~ 3 cV 2 . Allowing inde- 
pendent oscillations, we measure antineutrino parame 



|Ato 2 | = (2.50^;2s) x 10~ 3 cV 2 . A comparison of the 
neutrino and antineutrino mass splittings shows them 
to be in excellent agreement. These results provide the 
world's most precise measurement to date of these mass 
splittings for both neutrinos and antineutrinos. 
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